Multi-ion Mach-Zehnder interferometer with artificial nonlinear interactions 
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We show how to implement a Mach-Zehnder interferometry based upon a string of trapped ions 
with artificial nonlinear interactions. By adiabatically sweeping down/up the coupling strength 
between two involved internal states of the ions, we could achieve the beam splitting/recombination. 
Based on current techniques for manipulating trapped ions, we discuss the experimental feasibility 
of our scheme and analyze some undesired uncertainty under realistic experimental environment. 
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I. INTRODUCTION 

Due to long coherence time of some specific hyperfine 
states, high controllability of operations and high effi- 
ciency of state detection, the trapped atomic ion system 
has been considered as a promising candidate for quan- 
tum information processing. In particular, since they 
could be nearly perfectly prepared in some entangled 
states, the trapped ions have attracted considerable at- 
tention for high-precision quantum metrology [1-5]. 

It is generally believed that the measurement preci- 
sion can be enhanced from the standard quantum limit 
(SQL) or shot noise limit to the Heiscnbcrg limit by 
utilizing multipartite entanglement [6]. For example, 
it has demonstrated that the measurement precision 
may reach the Heiscnbcrg limit by using a NOON state 
(|7V) |0) b + |0) Q |7V) 6 ) /\/2, which is a superposition of 
N particles in mode a with zero particle in mode b and 
vice versa. However, most of the relevant proposals on 
high-precision interferometry of trapped ions are subject 
to limited numbers of ions or to the requirement for indi- 
vidual addressing of ions, which restricts the scalability 
of those interferometry schemes. 

In this article, we consider a Mach-Zehnder (MZ) in- 
terferometry using large numbers of entangled trapped 
ions. The MZ interferometry consists of a beam split- 
ter for splitting the incoming state and another beam 
splitter for state recombination. It has been shown that 
the adiabatic MZ interferometry is an optimal candidate 
for high-precision measurement [7, 8]. To achieve such 
an interferometry with trapped ions, we have to real- 
ize a nonlinear giant-spin Hamiltonian with anisotropic 
interaction. The key point is how to generate the nonlin- 
ear interactions and then achieve the beam splitting and 
recombination. By applying a pair of lasers to manip- 
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ulate independently their center-of-mass (COM) modes 
(cither transverse or longitudinal vibrational modes), the 
nonlinear giant-spin Hamiltonian with anisotropic inter- 
action could be simulated bases upon an array of trapped 
ions in a linear trap. By adjusting the coupling strength 
(i.e. the Rabi frequency) between the two involved hyper- 
fine states, the MZ interferometry could be carried out 
adiabatically. With current techniques for manipulating 
trapped ions, we discuss the experimental feasibility of 
our scheme and the difficulty due to decoherence in real- 
istic experimental environment. 

The paper is structured as follows. In the next section, 
we describe the physical system and derive its effective 
Hamiltonian. In Sec. Ill, we show how to realize the 
MZ interferometer by using quantum adiabatic processes 
through dynamical bifurcations. Then, we discuss the 
experimental feasibility in Sec. IV and analyze the un- 
desired uncertainty due to decoherence in Sec. V. In the 
last section, we briefly summarize our results. 



II. PHYSICAL SYSTEM AND ITS EFFECTIVE 
HAMILTONIAN 



We consider N identical ions confined in a linear trap, 
as shown in Fig. 1, where the three ionic states under 
our consideration are denoted by two hyperfine ground 
spin states (\\), and an excited state |e), respectively. 
In our scheme, all the ions are irradiated simultaneously 
by two traveling- wave laser beams with different frequen- 
cies, which coherently couple the hyperfine ground states 
through the optically stimulated Raman transitions by 
adiabatically eliminating the excited state |e). Here we 
assume the detuning of the optical fields from the elec- 
tronic resonance to be much larger than the excited state 
linewidth and the corresponding coupling strength. As 
indicated in Fig. 1(a), under an usual rotating wave 
approximation to the laser frequencies, the interaction 
Hamiltonian between the fields and the ions is written as 
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FIG. 1: (color online) Laser coupling scheme for simulating the Hamiltonian (2). (a) Physical realization of the nonlinear 
term of J^: a Stark shift regarding the level is created if two laser beams uniformly radiate the ions for the excitation 
from to |e)- with a large detuning A', (b) Physical realization of the terms of J x and J z : two additional laser beams are 
used to induce the carrier transitions corresponding to the term of J x ; and therefore the detuning to the resonant transition 
S = loo — (u' A — uj' b ) will generate the longitudinal-field term 8J Z . See more details in Appendix A. 
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where a\ = |4-)* (4-| — It)* (tl is the population inversion 
operator for the i-th ion, (a) is the creation (annihi- 
lation) operator of the center-of-mass (COM) mode and 
f2 2 is the uniform effective Rabi frequency, v is the vibra- 
tional frequency of the COM mode and A is the detuning 
from the vibrational frequency. r\ z denotes the Lamb- 
Dicke parameter fc z y '%/2Mv with the ion mass M and 
the laser wavevector difference k z = ks, z — &A,z- 

By tuning frequency differences and radiation direc- 
tions of the pair of lasers along with independent ma- 
nipulation of transverse or longitudinal COM modes, we 
may obtain following effective Hamiltonian under some 
approximations and canonical transformations, 



H, 



eff 



= 6J Z - B X J X - XJt, 



(2) 



where the collective angular momentum for all the ions 
are defined as Jk = Sili^fe/^ with k — (x,y,z), the 
longitudinal field 5 is the internal energy gap between 
the two spin states (|^), |t)), and the transverse field 
B x corresponds to the effective Rabi frequency, and the 
effective spin-spin nonlinear interaction is given as A = 
&Q 2 z r)1/(N A). More details of our derivation are given 
in Appendix A. These parameters 8, B x , and A could be 
controlled by Rabi frequencies, wavevector differences of 
the lasers and the detunings. 
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FIG. 2: (Color online) The energy spectra of the ground state 
|0) (dashed line) and the first excited state |1) (solid line) 
with respect to B x for different values of 5 and N. If <5 = 0, 
there is a bifurcation from degeneracy to non-degeneracy of 
the two lowest states |0) and |1) when B x increases, where 
(a) N=10 and (b) N=40. The beam splitting can be achieved 
by adiabatic passage through such a bifurcation. If S ^ 0, the 
degeneracy breaks down even in the weak limit of B x , where 
S = A/4, with N=10 (c) and N=40 (d). The other parameters 
are chosen as r) z = 0.1, A = 0.01 MHz, Q, z = 2ttx100 kHz 
(N = 10) and 2ttx200 kHz(iV = 40). 



III. MZ INTERFEROMETRY VIA ADIABATIC 
OPERATIONS 

The MZ interferometry consists of a beam splitter for 
splitting input states and another beam splitter for re- 
combining output states. Based upon the trapped-ion 
system (2), the two beam splitters for a MZ interferom- 
etry could be achieved by adiabatic quantum evolution 
through dynamical bifurcations. 



In the collective spin representation, the trapped-ion 
system is regarded as an ensemble of spin- 1/2 parti- 
cles. Let \jm) stand for the joint eigenstate of the SU(2) 
Casimir operators J 2 and J z , which satisfy the relations 

J 2 \ jm) = j(j + 1) \jm) , 
Jz \ jm) = m \jm) , 

with j = f and m = f ,f - 1, -f [10]. |f ) 
(|t> t)) means au i° ns m the hyperfine state \l) (|t))- 
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We first consider the case (5 = 0. In the strong coupling 
limit B x 3> |A|, the ground state |0) and the first excited 
state |1) for the Hamiltonian i? e // arc non-degenerated, 
where |0) = exp(i-| J y ) \ j = f , m = f) is a spin coherent 
state, and |1) is a superposition of different states \jm). 
When we adiabatically tune B x from the strong coupling 
limit (B x > |A|) to the weak coupling limit (B x < |A|), 
the ground state |0) and the first excited state |1) evolve 
into the states |f, — f ) an< ^ I"T't)' respectively. If 
B x = 0, the two lowest states turn to be degenerated. 
So in such an adiabatic passage, there exists a transition 
from non-degeneracy to degeneracy between |0) and |1), 
which can be regarded as a bifurcation [11, 12], as shown 
in (a) and (b) of Fig. 2. Similar to a single-particle 
MZ interferometry, where the ground state |0) and the 
first exited state |1) are utilized as two paths, the two 
states of maximum spin |f,— y) and |y,y) can a ^ so 
be used as two paths for an N-particle MZ interferome- 
try. Therefore, the achievement of the maximally path- 
entangled state |* p ) = 4= f ) + |f, -if)) implies 
accomplishment of the beam splitting. Then we turn off 
the lasers for a duration T so that the maximally path- 
entangled state \^f p ) evolves into 
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where ujq is the frequency difference between the two hy- 
pcrfine spin states \\) and 

There are two different methods for extracting the rel- 
ative phase information between the two paths of the 
N-particlc MZ interferometry. In the case of few trapped 
ions, this could be carried out using controlled-NOT 
gates C N OT^ est with the first ion being the control qubit 
and the rest ions being the target ones, followed by a 



Hadamard operation (Hi) on the first ion [3]. The rela- 
tive phase information can be obtained by measuring the 
population of the first ion in the state \\.) or [f). Tak- 
ing TV = 5 as an example, this individually addressing 
method can be summarized briefly as follows, 
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However, for systems of large numbers of ions, individ- 
ual addressing becomes challenging and it becomes very 
difficult to extract the relative phase with the controlled- 
NOT gate CNOT^ est . If we apply the inverse process of 
the first beam splitter for recombing the two paths, the 
relatives phase information indeed could be transferred 
into the population information of the two lowest states 
for the limit of B x ^> A. However, it is not easy to di- 
rectly distinguish these two lowest states. Fortunately, 
the two lowest states for symmetric (5 = 0) and asym- 
metric (8 ^= 0) systems of B x <C A are almost identical 
and the two lowest states (|f,— f), |f>f)) become 
non-degenerate for an asymmetric system (0 < S < A/4), 
see Fig. 2 (c) and (d). Therefore, to extract the relative 
phase, after applying the inverse process of the beam 
splitting process, we suddenly apply a proper nonzero 5 
(0 < 5 < A/4) and then adiabatically decrease B x from 
B x ^> A to B x = 0. This recombination procedure for ex- 
tracting the relative phase with two adiabatic processes 
could be summarized as, 



' \ adiabatic process of ( 5=0,B a= =0^-B a: »A COs(NuJ T/2) |0) — 1 Sm(NuJ T/2) |1) 

p ) ~ "~ * V2 
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Obviously, by measuring the population with all the ions 
in the state ||) or |t), the relative phase between the two 
paths could be obtained. 



IV. EXPERIMENTAL FEASIBILITY 

Our scheme is feasible with current ion trap tech- 
nology. If we employ a line of trapped 40 Ca + , the 
states ||), |t) and |e) are denoted by S 1 / 2 (to :) - = —1/2), 
Si/2(mj = 1/2), and Pi/ 2 (mj = -1/2) [13]. To adiabat- 



ically eliminate the ionic excited state |e), the detuning 
of the optical fields from electronic resonance should be 
much larger than the excited state linewidth and the cor- 
responding coupling strength, e.g. A' = 27rx20 GHz. On 
the other hand, generation of the nonlinear term A in 
the Hamiltonian (2) requires that the detuning A should 
be much smaller than the frequency v of the longitudi- 
nal COM mode. As a result, in our case we may take 
v = 1 MHz and A = O.OIia Besides, the coherence 
time of the motional states has reached 100 ms [14]. To 
change the effective coupling B x adiabatically, we may 
set B x (t) = a — f3t, where a and /3 are nonnegative pa- 
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ramcters, for example, a — 500 kHz and (3 = 50 kHz/ms 
with < t < 10 ms in the case of N = 40. Consider the 
system in a free evolution by T = 5 ms. After a duration 
for about 35 ms including three adiabatic processes plus 
the free evolution, we could obtain the relative phase by 
measuring the population of all the ions in state 1^' — t) 
or|f,+f>. 

Higher measurement precision for the trapped-ion MZ 
interferometer requires more ions involved. With more 
ions involved, however, besides more attentions should be 
paid to larger coupling strength, we also need to consider 
the detrimental influence from decoherence. Fortunately, 
currently available technologies enable the effective cou- 
pling strength to be on the order of MHz, which meets 
the requirement of our proposal. Therefore, we focus be- 
low on the influence from decoherence. 



V. ESTIMATION OF FREQUENCY 
UNCERTAINTY UNDER DECOHERENCE 



way is to enlarge A' by ten times, which leads to 7T to 
be about 0.001, and thereby the measurement precision 
would approach the Heiscnbcrg limit. Further improve- 
ment could be made by using refocusing pulses during 
the adiabatic operations to suppress the decay from |e). 
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There are several sources of decoherence. Here, as an 
example, we only consider the depasing decoherence in- 
fluence on the estimation of the ionic Zeeman splitting 
luq. The maximally entangled state of N trapped ions 
undertaking MZ interferometry decreases the measure- 
ment uncertainty of u>o by ^/N times compared to the 
trapped ions in product states [3, 5]. If decoherence due 
to environment is involved, the measurement precision 
would be seriously affected [15, 16] because entangled 
states are more sensitive to decoherence than product 
states. A recent experiment [17] has reported that the 
NOON state with 8 to 14 trapped ions subjects to cor- 
related dephasing noise mostly because of magnetic field 
fluctuation and shows super-decoherence, where the co- 
herence decays quadratically with the number of qubits. 
This correlated dephasing could be modeled by the Lind- 
blad equation in the free evolution, 

P = -iuo[J z , P\ + l&JzpJz - JzP ~ P J z) ( 6 ) 

where 7 is the correlated dephasing rate, which is A^ 2 7o 
with 70 the dephasing rate of a single qubit. Here, we 
consider 70 = Jq^I/A' 2 to be the effective two- level de- 
phasing rate [18] due to the two-photon Raman process, 
where 79 is the spontaneous decay rate of the excited 
state |e) for a single qubit. The uncertainty Aw/a>o 
shown in Fig. 3 is calculated in the equivalent spin rep- 
resentation, where the details are given in Appendix B. 

Repeating the implementation of our proposal by one 
hundred times, we find that the correlated dephasing seri- 
ously spoils the measurement precision of the entangled 
state, yielding the precision even lower than the SQL. 
By setting 7^ w 2tt x 20 MHz, il z = 2tt x 200 kHz, and 
A' = 2tt x 20 GHz, and considering free evolution by 
T = 5 ms in each implementation, we have 7T = 0.1 
and the uncertainty of uj in the case of three-ion entan- 
glement is larger than the SQL, see Fig. 3. In order to 
enhance the precision, we have to reduce 7. A possible 



FIG. 3: (color online) The frequency uncertainty Aw/wo verse 
the total number of ions N. The dotted curve represents the 
estimated uncertainty in SQL corresponding to the product 
state, and the dash-dotted curve is the Heisenberg limit (HL) 
with the NOON state. We take T = 5 ms and u « 2tt x 3 
GHz at the magnetic field B — 0.1 Tesla. 



VI. CONCLUSIONS 

In summary, we have shown how to adiabatically carry 
out a MZ interferometer based upon multiple trapped 
ions with artificial nonlinear interactions. Influence from 
decoherence has also been discussed in our treatment. 
Since our interferometry proposal involves large numbers 
of trapped ions and works by global operations, we argue 
that the proposal would be useful for high-precision quan- 
tum metrology with trapped ions toward the Heisenberg 
limit [19]. 

We should also emphasize that, although our discus- 
sion above has only focused on the case of A > and A > 
for a trapped-ion MZ interferometer, our model enables 
to study spin squeezing in an ion trap by setting A < 
for A < 0, which could be done by changing the directions 
of the laser beams k' A and k' B to make B x < [20]. Our 
model also enables to investigate interaction blockade of 
spin flip with trapped ions for weak B x [21]. 
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Appendix A: Generation of the effective 
Hamiltonian 



Below, we show how to generate the effective Hamilto- 
nian (2) in Section II. 

First, we discuss how to simulate the nonlinear term 
\J Z in the Hamiltonian H e jf- As indicated in Fig. 
1(a), irradiating the ions to the longitudinal COM mode 
by applying two non-copropagating laser beams with 
the wavevector difference k z along the z-axis, we have 
the Hamiltonian below in the rotating frame regarding 
Hq = (y — A)a'a under the rotating-wave approxima- 
tion, 
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Then applying a canonical transformation e s H2e with 
S = EL T7fH at " «H [ 22 ]. we obtain 



N 



(A2) 



Cooling the COM motion of the ions to its ground 
states, we obtain an effective spin-spin interaction in the 
Hamiltonian, 



(A3) 



with A = 8n 2 z r] 2 z /(NA) and J z = J^Zi °l/ 2 - 

The generation of the transverse- field term —B x J x re- 
quires additional lasers. As shown in Fig. 1(b), we set 
lj' a — uj' b = (Jo, implying a resonant Raman process with 
respect to the two ground spin states of the ions. Such 
operations yield the transverse-field term [22] 



see the schematic diagram Fig. 1(b). 

Since the manipulation on the transverse COM mode 
is independent of the longitudinal COM mode, the total 
effective Hamiltonian reads as 

H ef f =H 1 +H2+H 3 = SJ Z - B x J, - XJ Z . ( A6) 



Appendix B: Calculation of Aco 



For simplicity, we only consider the decoherence of 
the system during a free evolution. Suppose the sys- 
tem to be prepared in the maximally path-entangled 
statep p (0) = (|f,f) + |f,-f))«f,f| + (f,-f|)/2, 
where |^, ^} and |^, — ^) stand for all ions in |f) and 
respectively. After a time duration of T, the system 
evolves into 
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According to [16], the quantum Fisher information is 



given as 



F Q = Tr [p p {^,T)A 2 ] = N 2 T 2 e^ TN \ (B2) 

where the Hermitian operator A is the "symmetric 
logarithmic derivative" with its matrix elements de- 
^2 = -B x V" 4/2 = -B X J X , (A4) fined as Ay = 2[p' p (u) , T)] y /(pi + p s ), p' p (uJo, T) = 

— 'i— 1 f)n.-(t<)n T''\ I FfT wit/h rii nnH n- tViP picrpnva.lnps nf 



with B x an effective Rabi frequency [22, 23] and J x = 

EiIi4/2. 

If the Raman beams for generating the transverse-field 
term (—B X J X ) are non-resonant, i.e. the frequency dif- 
ference of the two lasers lo" a — ujg = loq — S with the 
detuning S, one can obtain the longitudinal-field term, 



dp p (u>o,T)/dT, with pi and pj the eigenvalues of 
Pp(uiO,T), conditional on pi + pj =0 and Aij = 0. 

The frequency uncertainty Alo satisfies the quantum 
Cramer-Rao bound, [24] 



Alu > 



Hx = 5 



al/2 = SJ Z 



(A5) 



in the measurements on a set of k probes. 
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